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SUMMARY 

Two  methods  for  estimating  the  three-parameters  of  the  Weibuil  extreme  value 
distribution  are  described  and  tested  using  Weibuil  distributed  data.  Graphical  techniques 
used  to  support  the  estimation  indicate  that  any  bias  may  arise  from  estimations  on  the 
boundary  of  the  parameter  space. 
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1.  INTRODUCTION 

^  No  adequate  computer  program  existed  at  ARL  to  estimate  the  parameters  of  a  three- 
parameter  Weibuil  distribution  from  a  set  of  experimental  data.  A  method  using  Maximum 
Likelihood  was  developed  and  programmed,  and  tested  with  Weibuil  distributed  data.  Contour 
plots  of  the  likelihood  function  have  been  made  as  a  check  on  the  convergence  of  the  solution, 
and  to  help  locate  the  maxima  when  the  method  fails  to  converge.  A  computerized  version  of  a 
graphical  method  using  extreme  value  probability  paper  [Ref  1  ]  was  also  developed,  and  the  results 
of  both  methods  compared.  , 


1.  THE  THREE-PARAMETER  WEIBULL  DISTRIBUTION 


The  three-parameter  Weibuil  distribution  is  an  extreme  value  distribution  which  takes  the 
following  form : 


F(x)  =  I  -exj 


and 


o) 

(2) 


with  the  conditions  v  >  *,  a  >  0,  F(t)  =  0,  F(v)  =  1—1/e.  Here  f(x)  is  the  probability  density 
function  of  the  random  variable  x,  a  the  dispersion  parameter,  v  the  characteristic  value  and 
*  the  lower  bound  of  x.  F(x)  is  the  cumulative  distribution  function : 


F(x)  = 


mix. 


3.  THE  MAXIMUM  LIKELIHOOD  METHOD 

The  likelihood  function  of  data  f(X{)  is  given  by 

eL',xt  =  flAXi) 

i- 1 

where  n  is  the  total  number  of  data  values  available.  Then,  from  (2) 

pl/ix) 


Taking  logarithms 
L 


W=n.n(-^)+£.nfcri-i(^y 

y>— «/  (-i  \v— »/ 

» /tin*— naln(v— »)+(a— 1)  J)  ln(Aj—  *)—  £  ( — — ) 
/- 1 


(3) 

(4) 
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If  we  make  the  substitution  V  =  v— «  and  reduce  the  magnitude  of  the  numbers  involved  by 
dividing  through  by  n,  we  get* 


L'(x\  1  * 

L(x)  - - =  lna— nln!/+(«  — 1)-  J)ln(JT|— *) — 

n  n  /  _  i 

The  Maximum  Likelihood  equations  for  V,  a  and  e  are: 


(5) 

(6) 

(7) 

(8) 


Equations  (6),  (7)  and  (8)  are  then  solved  for  V,  «.  and  «;  the  solution  method  employed  is  the 
second-order  Newton-Raphson  method  [Ref  3],  an  iterative  process  described  in  Appendix  1. 
The  method  requires  the  second  derivatives  to  provide  the  terms  for  the  Hession  matrix 
[32L/J0(Sfy].  These  are : 
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(9) 


To  check  that  a  maximum  is  being  obtained,  L  is  printed  after  each  iteration,  the  deter¬ 
minant  of  the  Hessian  matrix  is  found,  and  this,  as  well  as  the  terms  on  the  leading  diagonal, 
must  be  negative. 

To  enable  this  method  to  be  used  for  determining  the  parameters  of  a  two-parameter 
Weibull  distribution,  i.e.  «  =  0,  it  is  only  necessary  to  set  «  =  0  in  (6),  (7),  (8)  and  (9),  and  solve 
(A3)  with  i=y'  =  2  where  a i  =  F  and  02  =  «. 


4.  INITIAL  ESTIMATES  OF  THE  PARAMETERS  V,  «,  AND  « 


Initial  estimates  of  the  three  parameters  are  required  for  the  application  of  the  Newton- 
Raphson  method. 

From  equation  (7) 


1  n  In  /  Y. A  a  1  *  /  Yt A* 

-  I  ln(*-c)  =  lnv+l-  ln(*— 0— InF-  I [—-)  .  (10) 

«/-i  ni-i\  v  /  n/-i\  v  / 


*  Any  further  reference  to  “likelihood  function”  in  the  remainder  of  the  report  will  strictly 
mean  the  scaled  log  likelihood  function  L. 
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From  equation  (6) 


V  Fn(t,V  V  ) 

:  ,?M  ■ 

right-hand  side  of  (10)  gives 

*  »<-i\  v  )  "I- 1 

ion  of  the  three-parameter  Weibull  ] 

-ln(l-fU)))  = 


Substituting  (11)  into  the  right-hand  side  of  (10)  gives 

-  =  \  £  (^)W*-*)  -  -  £  !"(*<-«). 

*  »<-l\  v  )  »/-l 

Now,  from  (1),  the  definition  of  the  three-parameter  Weibull  probability  function  gives: 


Substitution  into  (12)  leads  to 


1  =  — -21n(*(— «)+-S(— ln(l  -F(Xi)) .  In( *,-«)).  (13) 

an  n 

Values  of  F(Xt)  can  be  approximated  from  the  data  by  listing  it  in  ascending  order  of 
magnitude  and  ascribing  probabilities  to  the  ordered  data  viz: 


1  -F,  ~ 


n+l-i 


This  approach  is  used  by  Gumbel  [Ref  3],  but  he  also  uses  a  correction  parameter  An,  to 
allow  for  different  sample  sizes  n,  to  obtain  a  better  estimate  of  1  —Ft.  The  corrected  value  is: 


where  An  =  for  n  <  100, 

nln(n-t-l)—  ln(«!) 

,  [ln(2wn)— 2] ,  _ 

«  1  +  —  for  n  >  100. 

On  using  this  corrected  value  (13)  becomes: 

-  -  l  £  ln(*— .)  -  —  £  ln/^f-'  W,-«).  (14) 

a  B(-i  n  i-i  \  rt+1  / 

The  initial  estimate  of  <  is  taken  as  the  data  minimum,  A"min,  i.e.  t  =  Xmm-  Using  this  value 
in  (14),  a  first  approximation  to  a  is  obtained. 

The  initial  estimate  of  V  is  obtained  from  the  mean  of  the  data,  i.e. 

V  -  -  £  X,-t.  (16) 

Hi- i 

[I  " 

-  £  (Xt—  •)*  I  ,an 

"i- 1  J 

initial  estimate  using  this  function  contains  uncertainties  inherent  in  the  initial  values  of  «  and  a. 
Using  the  population  mean  gives  quite  a  reasonable  estimate  of  V  and  it  is  independent  of  the 
other  parameters. 


5.  LEAST  SQUARES  OF  THE  REDUCED  VARIATE  y 
Define  a  new  variate 

y  -  -B[lD(jr-*)-ln(v-«)J- 

From  equation  (I), 

1  -/X*)  -  exp[— *-*J 
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or 

y  =  — ln(— ln(l  -0—F(x))). 

From  (17)  it  can  be  seen  that  y  and  hence  ln(— ln(l  0— F(x)))  is  a  linear  function  of  ln(x— *). 

In  general,  using  experimental  data,  only  the  x  values  are  available,  but  F(x)  may  be  approxi¬ 
mated  using  order  statistics  as  in  the  previous  section  (in  this  case  factor  A»  was  not  used). 


i.e. 


n+1  — i 

«+l 


Using  the  graphical  method  of  Reference  1,  yu  determined  using  ln(—  ln(l  —  F(Xi))),  is 
plotted  against  ln(Jf<— «),  and  by  adjusting  c,  a  least  squares  best  fit  straight  line  is  fitted  to  the 
data,  from  which  *  and  v  can  be  determined. 

To  obviate  the  need  to  estimate  c,  the  computational  procedure  solves  equation  (17)  for 
v,  a  and  <  using  a  nonlinear  least  squares  fit  of  against  Xt.  The  method  used  for  doing  this  is 
described  in  Appendix  1. 


6.  TEST  DATA 

Both  methods  were  tested  using  data  that  were  known  to  be  Weibull  distributed.  The  test 
data  were  produced  as  follows. 

From  equation  (1) 

-  loO-fW) 

and 

*  =  «+(«— «x— taO-fX*)))1'*. 

By  choosing  random  values  of  F(x)  between  0  and  1,  a  sample  of  x  can  be  obtained  that  will  be 
Weibull  distributed  with  parameters  v,  a  and  c.  In  this  way  samples  ranging  in  size  from  20  to 
1,000  were  obtained. 

The  random  numbers  were  generated  using  the  DEC  system  library  subroutine  RAN. 
RAN  is  not  dock-dependent;  it  will  reproduce  the  same  random  numbers  if  the  same  number  of 
calls  to  the  subroutine  are  made.  To  produce  different  sets  of  random  numbers,  a  chosen  number 
of  calls  to  the  subroutine  is  made  before  generating  the  data  set. 


7.  GRAPHICS  SUPPORT 

Both  methods  have  some  supporting  graphical  output,  to  help  indicate  the  success  of  the 
estimation. 

7.1  Graphical  Output  for  the  Maximus  likelihood  Method 

Contour  plots  of  the  likelihood  function  as  a  function  of  c  and  <*  can  be  produced,  together 
with  a  contour  plot  of  v  over  the  same  (*,a)  domain.  It  becomes  obvious  from  these  plots  where 
the  solution  is  located,  and  when  the  solution  method  fails  to  converge  they  can  provide  alterna¬ 
tive  etimates.  Figures  1  and  2  are  examples  of  the  contour  plots  produced. 

7.2  Graphical  oatpat  fir  leaat  aqaane  of  variate  y  amthod 

Since  this  method  determines  the  best  linear  fit  of  to  ln(2f<— «)  the  yt  is  plotted  against 
ln(Jlf«— «),  together  with  the  fitted  function.  An  example  is  shown  in  Figure  3. 


8.  RESULTS 

The  two  methods  were  applied  to  a  large  number  of  generated  samples.  One  set,  generated 
with  v  -  1  -0,  >  -  3-0  and  *  -  0-7,  demonstrates  the  effect  of  sample  size,  while  a  larger  set, 
with  u  -  10,  «  -  2- 5  and  «  -  0-7,  demonstrates  the  effectiveness  of  both  methods  when  only 
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a  small  sample  is  available.  Sample  sizes  of  1,000  were  also  generated  with  various  values  of  the 
three  Weibull  parameters. 

To  determine  the  effect  of  sample  size  on  estimates,  various  sizes,  from  20  to  1,000  were 
used,  with  five  different  samples  being  produced  for  each  size.  The  mean  values  of  each  parameter 
for  each  sample  size  are  shown  in  Table  1. 

Generally,  the  larger  the  sample,  the  greater  the  accuracy  of  estimation.  This  is  an  observa¬ 
tion  from  the  known  asymptotic  efficiency  of  maximum  likelihood  estimators,  and  it  probably 
indicates  that  the  least  squares  method  is  also  very  efficient.  There  appears  to  be  little  difference 
between  the  accuracy  of  the  two  methods  for  the  data  considered  in  Table  1. 

Extensive  testing  was  carried  out  with  a  typical  sample  size  of  20.  Thirty-six  different 
samples  were  produced  for  v  =  1-0,  «  =  2-5  and  c  =0-7.  (Each  sample  is  identified  by  a 
number  1ST  which  represents  the  number  of  calls  to  the  random  number  generator  prior  to 
producing  the  set  of  random  numbers.  This  enables  any  sample  to  be  reproduced.)  The  results 
are  shown  in  Table  2,  together  with  the  sample  identification  number  1ST  and  the  minimum 
value  in  each  sample,  A'min. 

The  maximum  likelihood  method  failed  to  converge  in  only  one  of  these  cases  (for  this  case, 
the  least  squares  method  returned  values  of  0-968,  1-152  and  0-790  for  v,  a  and  t  respectively). 

In  three  other  cases  both  methods  failed  to  produce  estimates  for  a  three-parameter  dis¬ 
tribution,  although  they  provided  estimates  for  a  two-parameter  distribution.  Table  3  shows 
three  different  data  samples  that  lead  to  different  types  of  solution  using  the  maximum  likelihood 
method. 

Case  32  had  a  successfully  converged  solution  with  v  =  0-96,  a  =  2-572  and  €  =  0-7076. 
The  contour  plots  of  the  likelihood  function  and  v  are  shown  in  Figures  1  and  2  respectively. 
The  location  of  the  maximum  in  the  likelihood  function  is  enclosed  within  the  0-95  contour, 
the  +  indicating  the  location  of  the  solution  values.  In  case  8,  the  Maximum  Likelihood  method 
failed  to  converge  to  a  solution.  The  contour  plot  of  Figure  4  shows  that  a  peak  does  not  exist 
within  the  *,«  domain.  The  likelihood  function  has  a  maximum  value  on  the  t  =  A'mm  =  0-804 
boundary  at  «  ~  0-8.  Since  no  converged  solution  exists,  the  best  estimates  of  v,  <*  and  r,  from 
Figure  4  and  the  corresponding  contour  plot  of  v,  Figure  5  are  0-94, 0-8  and  0-804  respectively. 

Case  9  would  only  converge  to  a  solution  if  t  was  set  to  zero  for  both  methods  of  estimation, 
i.e.  a  two-parameter  Weibull  distribution  appears  to  be  the  best  fitting  distribution  for  this  data. 
The  contour  plots.  Figure  6  and  7,  show  that  no  peak  exists  for  c  >  0,  and  that  the  maximum 
value  of  the  likelihood  function  lies  on  the  c  =  0  boundary  at  a  ~  9  ■  5  and  the  corresponding 
value  of  v  ~  1-069  (the  maximum  likelihood  method  returned  values  of  »  =  1-0693  and 
a  =  9-5259).  Considering  all  36  cases  of  Table  2,  the  following  means  and  standard  deviations 
of  the  three  parameters  are  obtained. 


Weibull 

parameter 

Values  used 
in  data 
generation 

Least  squares  of  variate  Y 

Maximum  likelihood 

Mean 

S.D. 

Mean 

S.D. 

V 

1  0000 

1-0068 

0-0300 

0  9958 

0  0305 

a 

2-5000 

4-3525 

5-4620 

3-2398 

2  9364 

€ 

0-7000 

0-5714 

0-2344 

0-6483 

0-2183 

The  mean  values  obtained  for  a  and  e  have  been  affected  greatly  by  the  three  cases  in  which 
t  was  assumed  to  be  zero,  cases  9,  26  and  35.  If  these  are  ignored,  the  mean  values  of  a  and  c 
are  much  closer  to  the  generating  parameters.  For  both  methods  there  is  a  large  scatter  in  the 
value  of  <*.  Ignoring  cases  9,  26  and  35  the  following  values  are  obtained. 


Weibull 

parameter 

Values  used 
in  data 
generation 

Least  squares  of  variate  Y 

Maximum  likelihood 

Mean 

S.D. 

Mean 

S.D. 

V 

1  0000 

1-005 

0  028 

0-993 

0-028 

a 

2-5000 

2-944 

1-653 

2-457 

1-129 

€ 

0-7000 

0-623 

0163 

0-707 

0-096 

The  Maximum  Likelihood  method  provides  better  estimates  for  a.  and  <;  both  methods  accurately 
estimate  v. 

When  all  the  sample  data  is  pooled,  providing  a  sample  size  of  720,  the  following  estimates 
for  v,  a  and  t  are  obtained  from  the  two  methods. 


Weibull 

parameter 

Values  used  in 
Data  generation 

Least  squares 
of  variate  Y 

Maximum 

likelihood 

V 

1-0000 

1-0029 

1-0018 

CL 

2-5000 

2-5761 

2-4800 

C 

0-7000 

0-6792 

0-6884 

Both  methods  were  tested  using  data  generated  with  various  other  values  for  the  three 
parameters,  and  the  results  are  shown  in  Table  4.  In  most  cases  the  Maximum  Likelihood  method 
gave  the  best  estimates  of  the  three  parameters. 


9.  DISCUSSION 

The  Maximum  Likelihood  procedure  requires  relatively  good  initial  estimates  of  *  and  a, 
otherwise  the  procedure  may  tend  to  diverge.  The  results  given  in  Table  2  reveal  that  the  mean 
value  of  e  is  10%  lower  than  the  mean  value  of  A'min  and  if  the  initial  estimate  of  c  is  slightly 
less  than  A'min,  the  procedure  is  more  likely  to  converge.  The  values  of  «  are  plotted  against 
A'min  in  Figure  8. 

It  was  also  found  that  introducing  a  relaxation  parameter  into  the  Newton-Raphson 
solution  procedure  improved  convergence. 

Of  particular  interest  are  the  contour  plots  of  the  likelihood  function.  All  have  a  relatively 
steep  gradient  along  the  «  =  A'min  boundary  with  a  “ridge”  running  from  *  =  A'min,  *  =  1-0 
to  a  point  on  the  c  =  0  boundary.  There  is  only  one  maximum  located  on  this  ridge. 

Although  the  contour  plots  give  the  impression  that  the  function  has  a  well-defined  “peak”, 
in  fact  the  “peak”  is  poorly  defined,  as  can  be  seen  from  the  values  of  the  contours;  the  ridge 
is  the  dominant  feature  in  the  e,a  domain. 

The  contour  plot  of  v  indicates  that  it  changes  slowly  throughout  the  «,«  domain  and 
explains  why  the  Maximum  Likelihood  procedure  is  relatively  insensitive  to  the  initial  value  of  v. 

Various  distributions  from  Table  2  with  parameters  determined  using  the  Maximum 
Likelihood  method  are  plotted  in  Figure  9.  Three  of  the  distribi  ‘:r>ns  are  for  the  data  sets  given 
in  Table  3.  Another,  case  17,  was  plotted  because  of  its  high  *  value  but  relatively  good  estimates 
of  v  and  «.  A  plot  of  the  parent  distribution  (y  =  1  0,  a  =  2-5  and  «  =0-7)  was  plotted  as  a 
comparison. 

The  interesting  feature  is  the  difference  between  the  distribution  of  case  32  and  the  parent 
distribution.  The  parameter  estimates  for  case  32  are  nearly  equal  to  the  parameter  of  the  parent 
distribution  yet  there  is  a  marked  difference  between  the  distributions. 
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10.  CONCLUSIONS 


The  three  parameters  of  a  Weibull  distribution  can  be  successfully  determined  using  a 
Maximum  Likelihood  approach.  The  method  can  be  used  with  as  little  as  20  data  values,  and  in 
the  event  that  the  method  fails  to  converge,  the  associated  graphics  enable  the  three  parameters 
to  be  estimated. 

When  a  solution  is  located  on  either  of  the  c  boundaries  of  the  parameter  space,  the  data 
used  is  probably  biased,  and  thus  caution  should  be  exercised  in  its  use. 

The  least  squares  of  the  reduced  variate  y  method,  developed  purely  as  a  comparison  with 
the  Maximum  Likelihood  method,  is  generally  less  accurate,  especially  when  a  is  large,  and  there 
may  be  inherent  errors  in  the  method  caused  by  the  need  to  assign  probabilities  to  the  data. 

Its  one  advantage  over  the  Maximum  Likelihood  method  is  that  the  initial  parameter  values 
do  not  have  to  be  as  accurately  determined  to  ensure  convergence. 
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APPENDIX  1 

The  Newtoa-Raphson  Method 

If  X  =  Xk  is  an  approximation  to  the  solution  X  =  A  of  f(x)  =  0,  then  the  sequence 

Xk+i  =  Xk  —  (A1) 

/  (*k) 

will  converge  quadratically  to  X  =  A  if  the  following  conditions  apply : 

(1)  for  Monotonic  convergence  F(Xa)f(X„)  >  0  and  f'(x),f’(x)  do  not  change  sign  in  the 
interval  (X0,A); 

(2)  for  Oscillatory  convergence  f(X0)f'(X0)  <  0  and  f'(x),f’(x)  do  not  change  sign  in  the 
interval  (X0,Xi),  X0  ^  A  ^  X\. 


If  8*  represents  the  change  in  Xk  after  each  iteration  then  (Al)  can  be  written  as 

f'(XK)&z  =  ~f(XK).  (A2) 

In  terms  of  the  Maximum  Likelihood  function  and  representing  V,  a  and  t  by  ai,  az  and  a3 
respectively,  sequence  (A2)  can  be  used  to  solve  equations  (6),  (7)  and  (8),  viz. : 


Z2L  „  ZL 

- oat  = - , 

ZatZaj  Zat 

or  expressed  in  matrix  form 


‘  Z2L 

Z2L 

Z2L 

Zai2 

ZaiZaz 

Z2L 

a2L 

Z2L 

d022 

dazZaz 

Z2L 

_  da&ai 

Z032 

j  =  1,3,  /  =  1,3  (A3) 


~ 

"  ZL  ' 

8ai 

5ai 

8a 2 

ZL 

=  — 

Zaz 

Saz 

ZL 

_  _ 

This  system  of  equations  is  solved  for  ai,  az  and  a3  using  a  subroutine  SOLEQU  contained 
in  the  ARL  DEC  10  library.  It  uses  Gaussion  elimination  with  the  pivot  selected  as  the  largest 
element  of  the  first  row  of  each  submatrix. 


I 


APPENDIX  2 

Nonlinear  Least  Squares  Fit  Method 

From  (17) 

yi  =  — ot(ln(A'l— «)— ln(u— t))+e(. 

The  residual  sum  of  squares  is  given  as 

£  =  £*2  =  £  b<+*(ln(A',-«)-ln(u-0))2.  ( 

i- 1  /- 1 

Suppose  we  have  trial  values  a0,  v0  and  e0.  We  shall  increment  these  by  An,  Av  and  At  respectively 
so  as  to  minimise  E  using  linear  approximations. 

For  simplification,  let  ft(Xt,aoj+Aaj)  =  — a(ln{;fr-f t)— ln(v— <))  where  aj  (j  -  1,2,3) 
replaces  v,  a  and  e  respectively,  a0j  replaces  the  initial  values  v0,  a0  and  t0  and  A aj  replaces  the 
incremental  values  Av,  Aa  and  At  then 

£  =  £ln-fi(Xi*4+AajP  <2 

i=  1 


t )E  n  d/i 

=  -2  Z 

d&aj  i  =  l  oOj  aQj,Xi 


neglecting  higher  order  terms. 


For  a  minimum  — —  =  0 


"  d/i ;  "  d/i 

•••  I  -  -  -yt  =  X  MXu  ac+Aa,) .  -J\  j  =  1 ,2,3 

i  =  ioaj  a0j.Xj  i=i  oaj 

These  equations  must  be  solved  for  A af  on  which  Zfdfifiaf  depends. 


MXt,a0,+Aat)  ~MX,,a0j)+  £  ^  .A a,. 

j  =  i  oa0j 


Substitution  into  (A®  gives 


i. ifi 

i  -  I  a0j 


■  yi  =  X  ¥(fi(X,,ao,)+  i  ^-Aa,) 

.A'/  i  =  iSa/\  j  „  i  ia0)  / 


transposing. 


V  „  f  tv  W  v  <>/«  v  <>f*  .  i,  "  ;>/<  A/t  . 

X  ~-(yt~fi(Xua0j))  =  I  ~I- —  •  Aa/  =  £  X  —  c~~‘ 

i  i  i>ai  i  - 1  So/  i  - 1  oa0]  j  =  i ;  - 1  id)  oa0f 


L 


In  matrix  form 


,-ii>ai5aoi  i-i  Jai  SOoz  i- 1  3a  1  Sa<>3 

Aoi 

«  D/< 

Zr^^i 

i=ioai 

^  A  £  A  A. 

i-i3a2  3aoi  f.i  Sa2Sa<>2  , -  ]  3a2  3aoS 

Aoj 

= 

"  2 >/< 

£^ [A.  £^<VL 

Ana 

"  2>/< 

i-i  SasSaoi  i_i  3aj  ia0>  i- 1  3as  Saoj  _ 

Z  (yi-fi) 

Li=idaa 

In  terms  of  equation  we  have 

^  =  ln(u— c)— Inf^i— t) 

ia 

dy  a 

<>U  U~— C 

"by  a  a 

St  (Xi—t)  (u—  f) 

and  the  matrix  equation  above,  in  terms  of  u,  a  and  c  becomes 


1=1 

£3yjy  £3y3j>' 
/=  l  Du Dot  Du  Dc 

Av 

»  3>-3y 
/-  i  3a  3u 

"  32y  »  3y3y 

i_i3a2  /=]3a3t 

A a 

= 

isMK'-))) 

£ 

;-i 3t  3u 

£  ^ 
i=i  3<3a  /«i3t2 

At 

S3MKS3)). 

This  system  of  linear  equations  is  solved  for  Av,  A«  and  At  and  the  initial  values  of  «,  a  and 
t  updated;  the  calculation  being  repeated  until  there  is  convergence  to  final  values  of  u,  aandt. 
The  initial  values  are  determined  as  described  in  section  4. 


TABLE  2 


Results  Obtained  for  Data  Samples  of  Size  20,  for  o  =  1  -0,  a  =  2  -5,  c  =  0-7 


1ST 

Least  squares  of  variate  y 

Maximum  likelihood 

No. 

V 

€ 

v  |  a 

€ 

1 

0 

0-7626 

0-9764 

3-5608 

0-5814 

0-9717  1  3-9643 

0-6096 

2 

20 

0-7811 

0-9646 

2-8545 

0-6952 

0-9561  2-1999 

0-7575 

3 

40 

0-7819 

0-9915 

1-6390 

0-7465 

0-9798  1-4563 

0-7744 

4 

60 

0-7400 

0-9651 

4-2318 

0-4762 

0-9712  5-2771 

0  4676 

5 

80 

0-7958 

1-0805 

2-2609 

0-6652 

1-0656  2-0486 

0-7381 

6 

100 

0-7720 

1 -0345 

3-3494 

0-5800 

1-0279  1  4-0224 

0-5934 

7 

120 

0-7572 

0-9872 

2-1532 

0-6683 

0-9754  !  1-9483 

0-7212 

8 

140 

0-8040 

0-9682 

1-5214 

0-7896 

0-9400t  0-8000 

0  8040 

9* 

160 

0-7574 

1-0775 

12-2768 

0-0 

1-0693  9-5259 

0-0 

10 

180 

0-7530 

1-0090 

6-3284 

0-3102 

10004  4-2105 

0-5825 

11 

200 

0-7241 

1-0006 

5-5321 

0-2806 

0-9909  1  3-9904 

0-5316 

12 

220 

0-7063 

1-0369 

4  •  5330 

0-3747 

1-0295  4-3204 

0  4904 

13 

240 

0-7866 

1-0065 

2-6827 

0-6893 

0  9967  2-1951 

0-7567 

14 

260 

0-8441 

1-0327 

1-5852 

0  8043 

1-0189  1-4248 

0-8349 

15 

280 

0  7612 

1-0089 

1-4850 

0-7233 

0  9934  i  1-4132 

0-7522 

16 

300 

0-7283 

0-9871 

1-6929 

0-6678 

0-9727  1-5384 

0  7154 

17 

320 

0-8056 

1-0264 

7-8281 

0-3614 

1-0196  4-5134 

0  6681 

18 

340 

0-7592 

0-9661 

1-2918 

0-7340 

0-9490  ;  1  2379 

0  7555 

19 

360 

0-8098 

1-0157 

1-6501 

0-7725 

1-0056  1-6338 

0  7965 

20 

380 

0-7661 

1-0062 

1-8961 

0-7003 

0-9895  1-6242 

0  7501 

21 

400 

0-7139 

1-0268 

2-1630 

0-6186 

1-0119  i  1-8577 

0  6856 

22 

420 

0-7274 

1-0507 

6-5505 

0-1565 

1-0359  3-5306 

0-5838 

23 

440 

0-8635 

1  -0330 

1-9446 

0-8197 

1-0235  !  1-7226 

0-8512 

24 

460 

0  8033 

1-0168 

1-8691 

0-7476 

1-0067  1-8798 

0-7774 

25 

480 

0-7233 

0-9987 

2-4818 

0  6186 

0  9865  2  0546 

0  6916 

26* 

500 

0  7541 

0-9853 

30-1486 

0  0 

09806  15-2848 

0  0 

27 

520 

0-7995 

1-0030 

2-6423 

0-7023 

0-9948  2-3492 

0  7576 

28 

540 

0-7539 

0-9740 

2-4077 

0  6705 

09652  2-1896 

0  7193 

29 

560 

0-8161 

1-0186 

4-6199 

0-5624 

1  0083  30710 

0  7343 

30 

580 

0  8325 

0-9816 

2-7091 

0-7309 

0-9744  2-3125 

0-7928 

31 

600 

0-8072 

1-0166 

1-8863 

0-7504 

1  0045  1  6600 

0  7910 

32 

620 

0-7772 

0-9715 

3-4026 

0-5923 

0  9600  2-5722 

0-7076 

33 

640 

0-7978 

1  -0419 

2-0191 

0-7278 

1-0044  1  9807 

0  7655 

34 

680 

0-7126 

1-0157 

3-2307 

0-5302 

1-0041  25730 

0  6479 

35* 

680 

0-7528 

1-0318 

17-1216 

0-0 

1-0249  107364 

0  0 

36 

700 

0-7452 

0-9631 

1-5082 

0-7119 

0  9512  !  1-5122 

0-7344 

*  Only  solutions  for  two-parameter  Weibull  distribution  existed, 
t  Solution  would  not  converge,  values  obtained  from  contour  plot. 

1ST — used  in  the  generation  of  random  numbers;  given  here  so  that  any  of  these  random 
number  data  sets  may  be  reconstituted. 


TABLE  3 


The  Three  Data  Sets  Used  in  Producing  the  Contour  Plots 
in  this  Report 


Case  32 

Case  8 

Case  9 

0-777221 

0-803982 

0-757371 

0-787285 

0-809143 

0-775983 

0-788506 

0-826658 

0-802988 

0-801296 

0-843383 

0-870959 

0-837540 

0-854666 

0-877282 

0-873750 

0-866742 

0  911713 

0-915978 

0-868289 

0-969554 

0-917450 

0-872045 

1-00788 

0-928643 

0-899308 

1-01243 

0-935790 

0-911857 

1-02508 

0-970270 

0-926751 

1-05588 

0-975081 

0-948090 

1-07683 

0-976064 

0-948586 

1  08666 

0-978354 

0-955043 

1-09314 

0-984281 

0-992642 

1-11991 

0-991846 

1 -05085 

1-12657 

1-00006 

1-05241 

1-13691 

1-00728 

1-11716 

1-16721 

1-03115 

1-22371 

1-17037 

1-15549 

1-23370 

1  20426 

Converged 

Did  not  converge 

Only  converged 

solution 

to  a 

solution 

to  a 

two-parameter 
solution, 
i.e.  «  =  0-0 

L 


TABLE  4 


'Estimates  of  the  Weibull  Parameters  Using  the  Two  Method  with  Various  Distributions  of  Sample 

Size  1000 


Weibull  parameters 
of  the  test  data 

Least  squares  of 
variety  solution 

Maximum  likelihood 
solution 

1ST 

V 

1 

«  1 

. 

V 

a 

t 

V 

-  !  * 

0 

100 

20 

90 

100  012 

80-354 

57-906 

100  003 

21  590  88-360 

2000 

100 

20 

90 

100  002 

68-970 

65-872 

99-997 

31-663  84-189 

0 

100 

20 

10 

100- 107 

80-360 

-278-869 

100-030 

21-613  -4-873 

0 

100 

5 

90 

100  037 

5-536 

88-390 

100  015 

4-832  89-597 

0 

100 

5 

10 

100-334 

5-536 

-4-492 

100-137 

4-832  1  6-371 

0 

l 

20 

0-9 

1-000 

80-356 

0-579 

1  000 

21-598  j  0-884 

0 

1 

3 

0-7 

1-001 

3-020 

0-682 

1-001 

2-842  I  0-694 

2000 

1 

3 

0-7 

1-000 

3-180 

0-683 

0-999 

3  096  j  0  691 

3000 

1 

3 

0-7 

1-009 

3-319 

0-683 

1-008 

3-181  1  0  693 

4000 

1 

3 

0-7 

0-999 

3-123 

0-689 

0-998 

3  084  j  0  694 

1000 

1 

3 

0-7 

1-004 

3-072 

0-701 

1-003 

2-785  0-717 

Dispersion  parameter,  a 


FIG.  1  CONTOUR  PLOT  OF 
(Data  are  shown  in  Ti 


Variate  V-,  =  —  Ln  (— Ln  (1.0  —  F(Xj))) 


FIG.  3  A  PLOT  OF  VARIATE  Y  FOR  A  DATA  SET  OF  500,  SHOWING 

THE  BEST  FIT  STRAIGHT  LINE  FITTED  BY  THE  LEAST  SQUARES 
OF  VARIATE  Y  METHOD 


Dispersion  parameter,  a 


I 

* 


Cut  off  e 


FIG.  5  CONTOUR  PLOT  OF  PARAMETER  u  FOR  CASE  8 
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jo  'jaiauiejed  uojsjadsjo 


Cut  off  e 


FIG.  6  CONTOUR  PLOT  OF  THE  LIKELIHOOD  FUNCTION  FOR  CASE  9 
(Data  ara  shown  in  Table  3) 


FIG.  8  MAXIMUM  LIKELIHOOD  ESTIMATE  OF  e  AGAINST  Xmin 
FOR  THE  36  CASES  SHOWN  IN  TABLE  2 


.1 


FIG.  9  WEIBULL  DISTRIBUTIONS,  f(X>,  DETERMINED  BY  THE 
MAXIMUM  LIKELIHOOD  METHOD  FOR  VARIOUS  DATA 
SAMPLES  CHOSEN  FROM  TABLE  2,  SAMPLE  SIZE  =  20 
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